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Abstract  A  novel  process  is  identified  whereby  equatorial  Rossby  (ER)  waves  maintain  warm  sea  surface 
temperature  (SST)  anomalies  against  cooling  by  processes  related  to  atmospheric  convection  in  the  western 
Indian  Ocean.  As  downwelling  ER  waves  enter  the  western  Indian  Ocean,  SST  anomalies  of  +0.1 5°C  develop 
near  60°E.  These  SST  anomalies  are  hypothesized  to  stimulate  convective  onset  of  the  Madden-Julian 
Oscillation.  The  upper  ocean  warming  that  manifests  in  response  to  downwelling  ER  waves  is  examined  in  a 
mixed  layer  heat  budget  using  observational  and  reanalysis  products,  respectively.  In  the  heat  budget, 
horizontal  advection  is  the  leading  contributor  to  warming,  in  part  due  to  an  equatorial  westward  jet  of 
80  cm  s-1  associated  with  downwelling  ER  waves.  When  anomalous  currents  associated  with  ER  waves  are 
removed  in  the  budget,  the  warm  intraseasonal  temperature  anomaly  in  the  western  Indian  Ocean  is 
eliminated  in  observations  and  reduced  by  55%  in  reanalysis. 


1.  Introduction 

The  temperature  of  the  oceanic  mixed  layer  (ML)  is  regulated  by  horizontal  advection,  radiative/heat  fluxes, 
vertical  advection,  entrainment,  mixing,  and  diffusion.  On  intraseasonal  time  scales,  the  tropical  ML  tempera¬ 
ture  (MLT)  varies  by  0.6°C  from  peak  to  peak  in  the  Indian  and  Pacific  Oceans  in  association  with  the 
Madden-Julian  Oscillation  (MJO)  [Madden  and  Julian,  1972;  Drushka  et  al.,  2012;  Halkides  et  al.,  2015]. 
Similarly,  intraseasonal  sea  surface  temperature  (SST)  variations  of  0.6°C  are  associated  with  the  MJO,  neglect¬ 
ing  diurnal  temperature  variations  that  can  exceed  2.0°C  [e.g.,  Mourn  et  al.,  2014;  de  Szoeke  et  al.,  2015; 
Ruppert  and  Johnson,  2014].  Surface  fluxes  are  primarily  responsible  for  patterns  of  SST/MLT  warming  and 
cooling  in  these  studies.  SST  is  observed  to  maximize  just  ahead  of  MJO  convection.  After  convection  begins, 
SST  rapidly  cools  and  reaches  a  minimum  ~5  days  later. 

However,  several  studies  have  observed  a  certain  class  of  MJO  events  that  deviate  from  the  previously 
observed  relationship  of  intraseasonal  convection  and  SST  anomalies.  For  these  select  MJO  events,  convec¬ 
tive  onset  is  observed  to  be  synchronized  with  downwelling  equatorial  Rossby  (ER)  waves  in  the  western 
Indian  Ocean  (10)  [Webber  et  al.,  2010,  2012b;  Rydbeck  and  Jensen,  201 7].  As  downwelling  ER  waves  enter 
the  western  10,  warm  SST  anomalies  phase  lock  with  enhanced  convection,  in  defiance  of  anomalous  cooling 
by  surface  fluxes.  Downwelling  ER  waves  and  the  associated  warm  SST  anomalies  are  hypothesized  to  initiate 
and  grow  MJO  convection  by  driving  convergence  in  the  atmospheric  boundary  layer  and  enhancing  surface 
heat  fluxes.  The  SST-driven  convergence  leads  deep  convection  by  several  days  and  is  hypothesized  to 
trigger  MJO  convection  [Rydbeck  and  Jensen,  201 7]. 

The  synchronization  of  intraseasonal  warm  SST  and  enhanced  convection  anomalies  during  MJO  initiation 
appears  to  be  unique  to  the  western  10.  Demott  et  al.  [2016]  observed  that  warm  SST  anomalies  in  the 
western  10  are  also  collocated  with  anomalously  increased  amounts  of  moist  static  energy.  Rydbeck  and 
Jensen  [2017]  showed  that  horizontal  convergence  of  moisture,  associated  with  downwelling  ER  waves,  is 
important  for  boundary  layer  moisture  buildup  that  precedes  MJO  initiation.  The  phase  relationship  of  warm 
SST  anomalies  and  enhanced  MJO  convection  evolves  as  the  MJO  propagates  eastward  in  the  10,  suggesting 
that  the  role  of  SST  anomalies  varies  as  a  function  of  the  MJO  lifecycle  [Gutzler  et  al.,  1994;  Sperber  et  al.,  1997; 
Webber  et  al.,  201 0,  201 2b;  de  Szoeke  et  al.,  2015;  Demott  et  al.,  201 6;  Rydbeck  and  Jensen,  2017]. 

In  this  study,  we  describe  a  process  by  which  warm  intraseasonal  SST  anomalies  develop  and  are  maintained 
in  the  western  10,  against  the  effects  of  cooling  by  surface  fluxes.  By  analyzing  the  ML  heat  budget,  we 
demonstrate  that  downwelling  ER  waves  induce  strong  horizontal  currents,  including  an  equatorial 
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westward  jet.  These  currents  drive  horizontal  advection  that  is  critical  to  the  growth  and  maintenance  of 
warm  SST  anomalies  prior  to  and  during  MJO  convective  onset.  Because  ER  waves  require  ~75  days  to  cross 
the  10  [Webber  et  at.,  2010,  2012a,  2012b;  Gottschalcket  at.,  2013;  Rydbeck  and  Jensen,  2017],  the  predictability 
of  MJO  events  associated  with  ER  waves  might  be  extended  by  a  similar  amount  of  time. 

2.  Data  and  Methods 

A  combination  of  observational  and  reanalysis  data  that  covers  a  20  year  period  from  January  1 993  to  December 
2012  is  used  to  examine  intraseasonal  SST  warming  associated  with  ER  waves.  Observed  sea  surface  height 
(SSH)  anomalies  are  obtained  from  Archiving,  Validation,  and  Interpretation  of  Satellite  Oceanographic  data 
(AVISO)  (0.25°)  [ Pascual  et  at.,  2006].  SST  from  the  National  Oceanic  and  Atmospheric  Administration  (NOAA) 
daily  optimum  interpolation  V2  that  includes  a  combination  of  advanced  very  high  resolution  radiometer 
(AVHRR)  and  Advanced  Microwave  Scanning  Radiometer-Earth  Observing  System  (AMSR-E)  are  used  in  obser¬ 
vational  composites  [Reynolds  et  at.,  2007].  The  AMSR-E/AVHRR  data  are  used  for  June  2002-October  201 1,  and 
the  AVHRR-only  data  are  used  for  the  remainder  of  the  1993-2012  period.  Currents  from  Ocean  Surface 
Current  Analysis  Real-time  (OSCAR)  (1/3°)  [Bonjean  and  Lagertoef,  2002,  Earth  Space  Research,  2009]  are  linearly 
interpolated  from  pentad  to  daily  data  for  use  in  compositing  with  other  daily  data.  OSCAR  currents  represent 
average  currents  in  the  upper  30  m  of  the  ocean.  TropFlux  surface  heat  and  radiative  fluxes  [Praveen  Kumar 
et  at.,  2012]  are  used  to  compute  the  net  fluxes  in  the  observed  ML  heat  budget.  NOAA  outgoing  longwave 
radiation  (OLR)  data  (2.5°)  [Liebmann  and  Smith,  1996]  are  used  as  a  proxy  for  atmospheric  convection. 

SSH,  temperature,  and  currents  from  the  Hybrid  Coordinate  Ocean  Model  (HYCOM)  Navy  Coupled  Ocean 
Data  Assimilation  (NCODA)  reanalysis  (1/12°)  [Metzger  et  at.,  2014]  are  also  used  for  comparison  with  observa¬ 
tions.  Climate  Forecast  System  Reanalysis  (CFSR)  [Saha  et  at.,  2010]  surface  fluxes  are  used  when  calculating 
the  ML  heat  budget  for  HYCOM  reanalysis.  ML  depth  (MLD)  is  computed  from  HYCOM  reanalysis  using  the 
fixed  density  criterion  of  0.125  kg  nrT3  referenced  from  the  surface.  The  MLD  is  used  in  the  calculations  of 
the  ML  heat  budget  in  both  HYCOM  and  observations.  These  results  are  not  significantly  different  when  other 
MLD  criteria  are  used. 

Downwelling  ER  waves  in  the  western  10  are  identified  using  the  method  of  Rydbeck  and  Jensen  [201 7].  To 
briefly  summarize  their  method,  SSH  anomalies  from  AVISO  are  filtered  in  wave  number-frequency  space 
to  remove  signals  outside  of  50-500  day  periods  and  planetary  zonal  wave  numbers  —1  to  —30  such  that 
ER  wave  signals  are  largely  isolated.  The  filtered  SSH  anomalies  are  averaged  from  6°S-6°N  and  59°E-69°E. 
The  resulting  time  series  is  used  as  an  ER  wave  index  for  compositing  periods  of  robust  downwelling  ER 
waves.  Events  greater  than  1.5  standard  deviations  are  used  for  compositing,  with  the  maximum  for  each 
event  set  as  lag  0  days.  The  dates  of  the  ER  wave  events  are  shown  in  Table  1  of  Rydbeck  and  Jensen 
[2017],  of  which  the  first  14  events  are  included  in  this  study.  Intraseasonal  anomalies  are  calculated  by 
removing  the  seasonal  cycle  and  applying  a  20-200  day  Lanczos  band-pass  filter  with  241  weights. 

3.  Equatorial  Rossby  Waves  and  Westward  Jets 

Composites  of  unfiltered  SSH  anomalies  and  currents  at  lags  —60  and  0  days  are  shown  (Figure  1 )  for  AVISO 
SSH  anomalies  and  OSCAR  currents  (left  column)  and  HYCOM  SSH  anomalies  and  surface  currents  (right 
column).  At  lag  —60  days  (Figures  la  and  1b),  an  eastward  jet  extends  across  much  of  the  10  in  both 
OSCAR  and  HYCOM  data.  The  eastward  jet  is  confined  along  the  equator  and  is  associated  with  an  eastward 
propagating  Kelvin  wave  that  is  downwelling  on  its  leading  edge  as  indicated  by  the  SSH  maxima  in  the 
eastern  10.  These  are  commonly  referred  to  as  Yoshida-Wyrtki  jets  and  are  forced  by  surface  westerly  winds 
[Yoshida,  1959;  Wyrtki,  1973].  Rydbeck  and  Jensen  [2017]  observed  that  the  downwelling  Kelvin  wave  in 
Figures  la  and  1b  is  forced  by  a  combination  of  the  semiannual  intermonsoon  and  MJO  westerlies.  In 
HYCOM,  the  magnitude  of  SSH  anomalies  associated  with  the  Kelvin  wave  in  the  central  and  eastern  10  is 
weaker  than  observations,  and  negative  SSH  anomalies  associated  with  upwelling  ER  waves  in  the  western 
10  are  stronger  than  observations.  This  is  partly  attributed  to  differences  in  the  climatology  used  to  calculate 
the  respective  SSH  anomaly  fields. 

After  reaching  Sumatra,  the  Kelvin  wave  is  partly  reflected  as  poleward  propagating  coastal  Kelvin  and 
westward  propagating  downwelling  ER  waves.  The  forcing  and  resonance  of  semiannual  Kelvin  and  ER  waves 
in  the  10  has  been  represented  in  numerous  models  of  varying  complexity  [i.e.,  Cane  and  Sarachik,  1981; 
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Figure  1.  Composites  of  unfiltered  SSH  anomalies  (shading;  cm)  and  currents  (vectors;  m  s_1)  are  shown  for  (a  and  b)  lag  -60  and  (c  and  d)  lag  0  days.  AVISO  SSH 
and  OSCAR  currents  are  shown  in  Figures  la  and  1c.  HYCOM  SSH  and  surface  currents  are  shown  in  Figures  1b  and  Id.  The  reference  vector  is  shown  in  the  bottom 
right  of  each  panel.  Vectors  with  magnitudes  less  than  0.05  m  s_1  are  omitted. 


Jensen,  1 993;  Han  et  al.,  1 999].  The  resultant  downwelling  ER  waves  are  associated  with  local  SSH  maxima  that 
flank  the  equator  near  4°N/S  and  arrive  in  the  western  10  at  lag  0  days  (Figures  1c  and  Id).  Positive  SSH 
anomalies  are  stronger  in  the  southern  hemisphere  than  the  northern  hemisphere  in  both  data.  However, 
HYCOM  reanalysis  slightly  underestimates  the  SSH  magnitude,  and  this  may  be  related  to  the  weaker 
downwelling  Kelvin  wave  at  lag  —60  days.  Along  the  equator  and  between  the  SSH  maxima,  an  equatorial 
westward  jet  is  present  at  lag  0  days.  Westward  jets  associated  with  downwelling  ER  waves  were  first 
observed  and  described  by  Delcroix  et  al.  [1991]  in  the  Pacific  Ocean.  Westward  jets  have  been  observed 
to  exceed  0.8  m  s_1  in  the  Indian  and  Pacific  Oceans  [ Delcroix  et  al.,  1991]  but  are  weaker  than  the  1  m  s_1 
and  greater  speeds  observed  for  Yoshida-Wyrtki  jets  in  the  10  [Wyrtki,  1973;  Knox,  1976;  McPhaden,  1982; 
Han  et  al.,  1999;  Mourn  et  al.,  2014;  Jensen  et  al.,  2015].  While  the  Yoshida-Wyrtki  jet  spans  most  of  the  10, 
the  zonal  extent  of  the  westward  jet  is  confined  to  the  half  wavelength  of  the  ER  wave. 

The  westward  transit  of  ER  waves  across  the  10  is  analyzed  using  composites  of  AVISO  and  HYCOM  SSH 
anomalies  averaged  from  6°S  to  6°N  (Figures  2a  and  2b).  Positive  SSH  anomalies  associated  with  downwelling 
ER  waves  originate  near  the  coast  of  Sumatra  and  propagate  westward  at  0.9  m  s_1,  consistent  with  the  first 
meridional,  first  baroclinic  mode  ER  wave.  After  weakening  in  the  central  10,  downwelling  ER  waves  reemerge 
in  the  western  10  with  increasingly  positive  SSH  anomalies.  The  arrival  of  downwelling  ER  waves  in  the 
western  10  is  coincident  with  the  local  development  of  negative  intraseasonal  OLR  anomalies,  as  previously 
observed  in  other  studies  [Webber  et  al.,  2010,  2012b;  Rydbeck  and  Jensen,  2017].  While  the  SSH  maxima  in 
HYCOM  are  -40%  weaker  than  observations  in  the  eastern  10  at  lag  —30  days,  the  maxima  are  comparable 
as  the  wave  arrives  in  the  western  10  at  lag  0  days. 

In  Figures  2c  and  2d,  a  westward  equatorial  jet  is  observed  to  propagate  along  the  path  of  downwelling  ER 
waves.  The  westward  jet  maximizes  in  the  central  and  western  10  between  lags  —25  and  0  days  and  is 
stronger  by  -30%  in  HYCOM,  in  part  due  to  HYCOM's  greater  resolution  [Joseph  etal.,  2012].  Near  the  leading 
edge  of  the  westward  jet,  between  lags  5  and  1 0  days  at  60°E,  HYCOM  MLT  and  NOAA  SST  anomalies  are  a 
local  maximum.  As  discussed  in  the  introduction,  this  local  SST  maximum  is  hypothesized  to  be  partly  respon¬ 
sible  for  MJO  convective  onset  by  driving  boundary  layer  convergence  and  enhancing  surface  heat  fluxes 
[Webber  et  al.,  2010,  2012b;  Rydbeck  and  Jensen,  2017].  The  observed  SST  maximum  is  0.1 5°C,  while  the 
corresponding  MLT  anomaly  in  HYCOM  is  0.1 2°C  (Figures  2c  and  2d).  The  HYCOM  SST  maximum  (not  shown) 
is  similar  in  magnitude  and  location  to  the  HYCOM  MLT  maximum,  suggesting  that  SST  may  be  used  as  a 
proxy  for  MLT.  This  relationship  is  utilized  to  compute  the  observed  ML  heat  budget  in  the  next  section.  In 
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Figure  2.  Composites  of  20-200  day  band-pass  filtered  (a)  NOAA  OLR  (shading;  W  m~2)  and  AVISO  SSH  (contours;  cm)  anomalies,  (b)  NOAA  OLR  and  HYCOM  SSH 
anomalies,  (c)  OCSCAR  zonal  current  (shading;  m  s_1)  and  NOAA  SST  (contours;  °C)  anomalies,  and  (d)  HYCOM  ML  zonal  currents  (shading;  m  s_1)  and  MLT 
(contours;  °C)  anomalies.  SSH,  SST,  and  OLR  anomalies  are  averaged  from  6°S  to  6°N.  Zonal  current  anomalies  are  averaged  from  3°S  to  3°N.  Phase  lines  of 
0.9  m  s_1  westward  and  5.0  m  s_1  eastward  are  indicated  by  bold  black  lines.  Positive  (negative)  contours  are  solid  (dashed),  and  the  contour  interval  is  shown 
at  the  bottom  right  of  each  panel.  The  zero  contour  is  bold. 


observations,  the  SST  maximum  and  OLR  minimum  near  60°E  are  largely  in  phase.  Similar  phasing  of  the  SST 
maximum  and  OLR  minimum  is  also  observed  in  HYCOM.  The  continued  warming  of  SST  in  the  western  10 
after  convective  onset  in  both  data  is  a  remarkably  interesting  feature  of  the  ER  wave  composites.  To 
further  examine  the  process  by  which  intraseasonal  SST  anomalies  warm  in  the  presence  of  enhanced 
MJO  convection,  composite  ML  heat  budgets  for  observations  and  reanalysis  are  analyzed,  respectively. 


4.  Mixed  Layer  Heat  Budget 

The  hypothesis  by  which  robust  downwelling  ER  waves  trigger  MJO  convective  onset  hinges  on  the 
emergence  of  warm  SST  anomalies  in  the  western  10.  These  anomalies  are  critical  to  the  SST  configuration 
that  accounts  for  up  to  45%  of  the  intraseasonal  atmospheric  boundary  layer  convergence  [Rydbeck  and 
Jensen,  2017].  To  determine  whether  SST  anomalies  in  the  region  are  generated  as  a  result  of  horizontal 
advection,  surface  fluxes,  processes  related  to  the  vertical  velocity  such  as  vertical  mixing  or  Ekman  pumping, 
or  any  combination  thereof,  an  ML  heat  budget  is  calculated  using  observations  and  reanalysis,  respectively, 
to  better  approximate  the  range  of  possible  solutions.  Because  observations  of  the  ML  are  scant  over  the 
region  and  period  of  interest,  ML  currents  and  MLT  are  estimated  using  OSCAR  currents  and  NOAA  SST.  In 
HYCOM  reanalysis,  average  ML  current  and  temperature  are  utilized. 
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The  intraseasonal  ML  heat  budget  is  written  as 

(^r)  =  -((7-VMLT)'  +  (Qs)'  +  R  (1) 

In  this  equation,  primed  terms  are  band-pass  filtered  for  20-200  days  to  capture  the  intraseasonal  signals 
appropriate  for  ER  wave  and  MJO  time  scales.  The  time  tendency  of  the  MLT,  (  a“LT)  ,  is  determined  by 
horizontal  advection,  -((7- VMLT),  scaled  net  surface  flux,  (Qs) ,  and  the  residual,  R.  The  net  surface  flux 
includes  latent  heat,  sensible  heat,  longwave  radiative,  and  net  shortwave  radiative  fluxes  and  is  scaled 
by  the  MLD  and  volumetric  heat  capacity  of  seawater.  Net  shortwave  radiative  flux  is  determined  as 
the  incoming  shortwave  radiative  flux  less  the  portion  that  penetrates  the  base  of  the  ML,  as  calculated 
in  Pacanowski  and  Griffies  [1998].  Confidently  calculating  processes  related  to  vertical  motions  is  difficult 
using  reanalysis  and  not  achievable  using  observations  for  the  region  and  time  periods  of  interest.  As  a 
result,  the  residual,  R,  includes  temperature  tendencies  due  to  Ekman  pumping,  entrainment,  vertical 
mixing/diffusion,  and  any  errors  in  the  calculation  of  the  other  terms.  The  ML  heat  budget  terms  are 
averaged  from  6°S-6°N,  58°E-62°E  (see  box  in  Figure  4)  to  capture  the  warming  coincident  with  MJO 
convective  initiation. 

In  the  observations,  the  MLT  increases  by  0.28°C,  from  —  0.1 6°C  at  lag  —23  days  to  0.1 2°C  at  lag  +8  days 
(Figure  3a).  Similar  to  the  observations,  the  MLT  in  HYCOM  increases  by  0.25°C,  from  — 0.14°C  at  lag  —23  days 
to  0.1 1  °C  at  lag  +5  days.  In  both,  the  warm  MLT  anomaly  is  sustained  for  several  weeks.  The  largest  contribu¬ 
tor  to  MLT  warming  in  both  budgets  is  horizontal  advection,  with  a  maximum  value  of  0.012  K  d_1  at  lag 
+1  day  in  observations  and  0.011  K  d_1  at  lag  —3  days  in  HYCOM.  The  respective  contributions  by  zonal 
and  meridional  advections  are  comparable  (not  shown).  Net  surface  fluxes  from  TropFlux  and  CFSR  contri¬ 
bute  to  MLT  warming  prior  to  convection  but  act  to  cool  MLT  after  convection  begins  near  lag  +2  days 
(see  Figures  2a  and  2b).  Between  lag  —10  and  +20  days,  the  residual  in  observations  mostly  contributes  to 
MLT  cooling.  The  residual  in  HYCOM  also  contributes  to  cooling  between  lag  +3  and  +15  days.  However, 
the  residuals  in  observations  and  HYCOM  are  out  of  phase  prior  to  lag  +2  days,  suggesting  different  roles 
for  unresolved  processes  at  the  leading  edge  of  downwelling  ER  waves. 

The  prominent  role  of  horizontal  advection  in  both  budgets  suggests  that  circulations  driven  by  ER  waves 
are  responsible  for  much  of  the  intraseasonal  warming  in  the  western  10.  To  quantify  the  magnitude  of 
the  warming  resulting  from  advection  by  intraseasonal  currents,  of  which  ER  waves  are  the  leading 
contributor,  the  horizontal  advection  term  is  separated  into  contributions  from  intraseasonal  and 
nonintraseasonal  currents: 


—  (L/-VMLT)'  =  — (V-VMLT)  -(Ju  -  l/)-VMLT)  (2) 

Primed  quantities  indicate  intraseasonal  anomalies.  The  term  on  the  left-hand  side  of  equation  (2), 
-((7- VMLT),  is  the  horizontal  advection.  The  first  term  on  the  right-hand  side  (RHS)  of 
equation  (2),  —((7  ■  VMLT)  ,  is  the  horizontal  advection  resulting  from  the  intraseasonal  currents.  The  second 
term  on  the  RHS  of  equation  (2),  —(((7—  (7)  -  VMLT)  ,  is  the  horizontal  advection  by  nonintraseasonal  currents. 
Nonintraseasonal  currents  include  high-frequency  variability  less  than  20  days  and  low-frequency  variability 
greater  than  200  days. 

Figures  3c  and  3d  show  select  terms  of  the  ML  heat  budget  with  (dashed)  and  without  (solid)  horizontal 
advection  by  intraseasonal  currents,  —((7  -  VMLT)  .  In  observations,  horizontal  advection  and  the  MLT 
tendency  are  weakened  by  -70%  at  lag  0  days  when  intraseasonal  currents  are  removed  (Figure  3c). 
The  contributions  to  horizontal  advection  by  the  zonal  and  meridional  intraseasonal  currents  are  compar¬ 
able  (not  shown),  indicating  that  the  meridional  currents  are  as  important  as  the  westward  jet.  The  dra¬ 
matic  reduction  of  the  MLT  tendency  results  in  an  MLT  maximum  that  is  weakened  to  near  zero  and 
shifted  to  earlier  lags  by  3  days.  From  these  results,  the  role  of  intraseasonal  currents  associated  with 
ER  waves  is  to  warm  the  MLT  and  shift  it  to  later  lags.  This  method  does  not  remove  all  effects  of  ER 
waves  on  the  MLT  tendency  and  anomaly.  For  example,  processes  such  as  vertical  advection  and  mixing 
by  ER  waves  remain. 
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Figure  3.  Intraseasonal  ML  heat  budget  terms  averaged  from  6°S-6°N,  58°-68°E  are  shown  for  (a  and  c)  observations  and 
(b  and  d)  HYCOM  reanalysis.  In  Figures  3a  and  3b,  all  terms  in  equation  (1)  are  shown.  In  Figures  3c  and  3d,  horizontal 
advection  by  intraseasonal  currents  is  removed  from  the  MLT  tendency,  horizontal  advection,  and  MLT  and  is  shown  in 
solid  lines.  The  dashed  lines  indicate  the  corresponding  complete  terms  shown  in  Figures  3a  and  3b.  In  each  panel,  MLT 
anomalies  are  multiplied  by  10_1. 


When  intraseasonal  currents  are  removed  in  HYCOM  (Figure  3d),  the  maximum  of  horizontal  advection  and 
MLT  tendency  are  both  reduced  by  -40%.  Without  horizontal  advection  by  intraseasonal  currents,  the  MLT 
maximum  is  reduced  by  55%  and  occurs  2  days  earlier.  As  in  observations,  the  results  from  HYCOM  suggest 
that  intraseasonal  ML  currents  associated  with  downwelling  ER  waves  are  critical  to  the  amplification  of  warm 
SST  anomalies.  Seiki  et  at.  [2013]  observed  during  the  CINDY/DYNAMO  field  campaign  in  the  10  that  horizon¬ 
tal  advection  by  oceanic  ER  waves  was  important  to  the  generation  of  intraseasonal  SST  variability  responsi¬ 
ble  for  modulating  MJO  initiation.  The  ML  heat  budget  of  Seiki  et  al.  [201 3]  revealed  that  net  surface  fluxes 
and  vertical  entrainment  could  not  fully  account  for  the  observed  SST  variability,  similar  to  the  results  shown 
in  Figure  3. 

Figure  4  shows  the  20-200  day  band-pass  filtered  current  anomalies  and  total  SST/MLT  from  observations 
and  HYCOM,  respectively,  at  lag  0  days.  The  currents  shown  in  Figure  4  are  those  that  are  suppressed  in 
Figures  3c  and  3d.  In  both  composites,  anomalous  westward  currents  with  maxima  of  0.3  m  s_1  are  present 
across  most  of  the  equatorial  western  10  and  predominantly  located  south  of  the  equator.  The  intraseasonal 
currents  associated  with  downwelling  ER  waves  act  on  the  temperature  gradient  to  yield  warm  anomalies  in 
the  western  10,  even  in  the  presence  of  cooling  by  surface  fluxes.  The  SST  gradient  is  unique  to  the  western  10 
and  likely  explains  the  preference  for  warm  SST  anomalies  associated  with  ER  waves  to  emerge  there  (see 
Figures  2c  and  2d). 
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Figure  4.  Composites  of  (a)  20-200  day  band-pass  filtered  OSCAR  current  anomalies  (vectors;  ms-1)  and  unfiltered  NOAA 
SST  (shading;  °C)  and  (b)  20-200  day  band-pass  filtered  HYCOM  ML  current  anomalies  (vectors;  m  s_1)  and  unfiltered 
MLT  (shading;  °C)  are  shown  for  lag  0  days.  Reference  vectors  are  shown  in  the  bottom  right  of  each  panel.  The  averaging 
box  for  the  budget  calculations  is  marked  in  white.  Vectors  less  than  0.025  m  s_1  are  omitted. 

5.  Discussion 

Processes  regulating  intraseasonal  SST  variability  in  the  western  10  are  examined  during  periods  of  strong 
downwelling  ER  waves.  SST  warming  associated  with  the  arrival  of  downwelling  ER  waves  in  the  western 
10  is  hypothesized  to  be  responsible  for  MJO  convective  onset  by  modulating  SST  gradients  that,  in  turn, 
generate  convergence  in  the  atmospheric  boundary  layer  and  enhance  surface  heat  fluxes  [Webber  et  al., 
2010,  2012b;  Rydbeck  and  Jensen,  2017].  Within  this  framework,  identifying  the  processes  responsible  for 
SST  warming  enhances  our  understanding  of  MJO  convective  onset  and  may  extend  current  MJO 
predictive  capabilities. 

Downwelling  ER  waves  in  the  10  are  linked  to  westward  propagating,  positive  SSH  anomalies  located  at  4°N/S 
and  a  westward  jet  along  the  equator.  As  the  dowelling  ER  waves  arrive  in  the  western  10,  warm  SST  anoma¬ 
lies  of  0.1 2-0.1 5°C  emerge  near  the  leading  edge  of  the  westward  jet.  The  phase  relationship  between 
convection  and  SST  anomalies  during  MJO  initiation  in  this  study  is  unique.  Warm  SST  anomalies  are  predo¬ 
minantly  in  phase  with  enhanced  convection  during  convective  onset  in  the  western  10.  Because  the  net 
surface  flux  is  negative  once  convection  begins,  horizontal  advection  is  largely  responsible  for  the  continued 
growth  and  maintenance  of  warm  SST  during  this  period. 

The  results  of  ML  heat  budgets  for  observations  and  reanalysis  suggest  that  continued  SST  warming  after 
convective  onset  is,  to  a  large  degree,  accomplished  by  horizontal  advection  driven  by  currents  associated 
with  downwelling  ER  waves.  These  results  support  the  remote  forcing  of  warm  SST  anomalies  in  the  western 
10  by  downwelling  ER  waves  and  emphasize  the  need  for  properly  resolving  ER  wave  propagation  and 
reemergence  in  coupled  models.  Horizontal  advection  derived  from  OSCAR  currents  should  be  interpreted 
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with  caution  as  the  meridional  currents  are  in  poorer  agreement  than  zonal  currents  when  compared  to 
observations  [i.e.,  Johnson  et  at.,  2007;  Sikhakolli  et  a!.,  2013].  This  may  partly  account  for  discrepancies  in 
the  horizontal  advection  and  residual  terms  of  the  respective  budgets. 

Robust  downwelling  ER  waves  in  the  western  10  are  linked  to  the  generation  of  downwelling  Kelvin  waves  in 
the  western  10  75-90  days  earlier  [Webber  et  al.,  2010, 2012b;  Rydbeck  and  Jensen,  2017],  The  time  required  for 
ER  waves  to  transit  the  10  suggests  that  MJO  events  initiated  by  downwelling  ER  waves  are  predictable  at 
10-13  week  leads.  Such  lead  times  greatly  exceed  the  current  3-6  week  forecast  capability  of  the  MJO 
[e.g.,  Neena  et  al.,  2014;  Kim  et  al.,  2014].  These  results  emphasize  the  need  for  fully  coupled  atmosphere- 
ocean  modeling  systems  to  gain  predictive  skill  of  MJO  initiation.  While  not  all  MJO  events  that  initiate  in 
the  western  10  are  associated  with  downwelling  ER  waves,  we  suppose  that  any  process  capable  of  growing 
and  maintaining  warm  SST  anomalies  against  cooling  by  the  net  surface  flux  might  also  stimulate 
MJO  convection. 
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